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Abstract — The MIMO wireless channel offers a rich ground for 
quality of service analysis. In this work, we present a stochastic 
network calculus analysis of a MIMO system, operating in 
spatial multiplexing mode, using moment generating functions 
(MGF). We quantify the spatial multiplexing gain, achieved 
through multiple antennas, for flow level quality of service (QoS) 
performance. Specifically we use Gilbert-Elliot model to describe 
individual spatial paths between the antenna pairs and model the 
whole channel by an iV-State Markov Chain, where N depends 
upon the degrees of freedom available in the MIMO system. 
We derive probabilistic delay bounds for the system and show 
the impact of increasing the number of antennas on the delay 
bounds under various conditions, such as channel burstiness, 
signal strength and fading speed. Further we present results for 
multi-hop scenarios under statistical independence. 

I. Introduction 

The last decade has seen an explosive growth in wireless 
network services and consequently an increasing demand of 
high data rates on wireless networks. Multiple-input-multiple- 
output (MIMO) |[T], |2j communication has generated a lot 
of interest in recent years for providing such high data rates. 
It offers enormous capacity gains as compared to traditional 
single input single output (SISO) communication J3J. The 
importance of MIMO technology can be seen from the fact 
that IEEE 802.11 is actively pursuing 802. lln which promises 
high throughput using multiple antennas. To date, a significant 
amount of work on MIMO has been done in quantifying its 
capacity limits under different channel conditions, antenna 
settings, etc., using classical information theory under the 
assumption of infinite delay [4|. This assumption has been 
working well for delay insensitive applications and for chan- 
nels without memory. 

However, recently there has been a rapid growth of delay 
sensitive applications such as VoIP and IP video, which are 
also run on wireless networks where MIMO is being deployed. 
In addition, MIMO links have been considered for delay 
critical applications such as MANETS [5|, for they overcome 
the problems of low system throughput. Hence, it is also 
essential to provide quality of service (QoS) (such as delajQ 
bounds) on MIMO links in addition to high data rates. For 
this, a critical issue is to analyze their flow level performance. 

The focus of this paper is on flow-level delay analysis of 
a spatial multiplexing MIMO wireless channel. Specifically, 
we derive probabilistic delay bounds for the channel and 
show the impact of increasing the number of antennas on 

'Throughout the paper the term delay refers to queuing delay. 



the delay bounds under various conditions, such as channel 
burstiness, signal strength and fading speed. The theory, which 
the analysis is based on, is moment generating function (MGF) 
based stochastic network calculus. In stochastic network cal- 
culus, a communication system is described by the arrivals at 
the ingress, the service it provides and the departures at the 
egress. The stochastic network calculus framework provides 
the freedom of a general separation in the characterization of 
arrivals and service, which allows the combination of different 
arrival and service models. It enables the derivation of flow 
level performance bounds such as Prob [Delay > d] < e, 
i.e. the probability of the delay exceeding a certain threshold 
d is upper bounded by e. It also incorporates the ability to 
characterize multiplexing and demultiplexing of flows and the 
derivation of end-to-end performance bounds through compo- 
sition |6|, [7|. The MGF based stochastic network calculus, as 
proposed in |8 1, established and developed in [7|, enables more 
easily the consideration of statistical independence between 
arrivals and system service and between the service provided at 
consecutive hops. A comprehensive view of stochastic network 
calculus can be found in (9j, (Toj. 

In the literature, many works, which are related to QoS 
analysis of wireless channels, can be found. They include link- 
layer channel modeling JTT), where the focus has been on 
developing models for channels with memory, which capture 
the bursty nature of errors on wireless channels. Effective 
bandwidth theory has been employed to analyze wireless 
channels [12|. The concept of effective capacity, which is 



directly related to effective bandwidth, was developed in 1 13 1 



for wireless channel QoS analysis. By using large deviations 
techniques a similar approach was developed in fl4[ to relate 
effective bandwidth theory to channel capacity and derive de- 
lay bounds in conjunction with Markov modeling. In addition, 
stochastic network calculus, which as a more general theory 
covers effective bandwidth and effective capacity, has recently 
been applied for QoS analysis of wireless channels fT3)-p7). 
However, none of these works analyzed MIMO channels. 

There are two transmission techniques for MIMO systems. 
Spatial multiplexing is used to increase throughput by trans- 
mitting independent data flows over different antennas. Spatial 
diversity, on the other hand, increases reliability by transmit- 
ting redundant flows of information in parallel. There has been 
some works reported for MIMO diversity systems (18), fl9) . 
In p"8) , the concept of effective capacity has been used to 
carry out QoS analysis of MIMO diversity systems using finite 



state Markov chains to describe the received signal to noise 
ratio (SNR). However, to the best of our knowledge, QoS 
analysis of MIMO systems in spatial multiplexing mode is 
still missing. In addition, it is known that, besides average 
received SNR, the MIMO channel matrix, which describes 
the scattering environment, has a profound effect on channel 
capacity (3) and hence the delay. 

The contribution of this work is as follows. (1) We present 
QoS analysis of a MIMO system under spatial multiplexing 
mode. In the analysis, we incorporate statistical independence 
of arriving traffic and the service provided by the channel. We 
also present results for end-to-end performance in a multi-hop 
scenario, where statistical independence of consecutive hops 
holds. (2) The effect of individual channel matrix elements on 
delay bounds is studied keeping the received SNR constant. 
(3) We provide a method of aggregating the states based on 
degrees of freedom available in the MIMO system, which 
results in reduction of the state space from 2 L to L + 1. 

The paper is organized as follows. We first present the 
system model and the assumptions in Sec.[llJ and the necessary 
background on stochastic network calculus and MIMO in 



Sec. II-A and Sec. II-B respectively. We then carry out the 



analysis in Sec. Ill Numerical results for a parameterization 



according to IEEE 802. lln are presented in Sec. |IV| 

II. System Model, Assumptions and Preliminaries 

We consider the single user MIMO case of a transmitter and 
a receiver with N and M antennas respectively. The MIMO 
system adopts spatial multiplexing. We assume a discrete time 
block fading model in which the channel is constant over the 
time length of one data block. We introduce a service model 
based on Gilbert-Elliott (GE) channels (30), (21]. Each fading 
channel is modeled by two states, {good} and {bad}, where a 
good channel is assumed for an error free transmission while 
a bad channel represents an unsuccessful transmission. We use 
MGF stochastic network calculus (NetCal) to carry out QoS 
analysis of the above mentioned MIMO system. 

A. MGF Stochastic Network Calculus 

Stochastic network calculus is a framework for QoS analysis 
and flow level performance evaluation of communication sys- 
tems. This framework has been firstly used for characterizing 
fading channels in (15). We resort to a stochastic NetCal 
approach based on moment generating functions, introduced 
in (7), for the analysis of the considered MIMO wireless fading 
channel. The first analysis of block fading channels using the 
MGF based stochastic NetCal was carried out in |16|. We 
carry forward the properties of the stochastic network calculus 
approach to MIMO wireless channels. 

The stochastic NetCal with MGFs allows efficiently exploit- 
ing statistical independence of arrivals and service processes 
and the statistical independence of service provided at con- 
catenated hops. It enables the analysis of a wide variety of 
scenarios for different arrival and service processes where the 
MGF can be calculated. 



We build on the dynamic server definition from [8] to model 
the time varying capacity of the channel. The discrete time ar- 
rivals A(0, t) and service S(s, t) are assumed to be statistically 
independent and stationary random processes. S(s,t) is the 
service and A(s, t) denote the cumulative arrival data provided 
in the interval (s, t] respectively. The MGF of a stationary 
process X(t) is given by 



M 



,t) = E 



,8X(t) 



(1) 



while M x (0,t) = M x (-6,t) for parameter 9 > 0, t > 
and E [•] denotes the expectation. For a set of useful MGF 
properties we refer the interested reader to (7), (22). Using 
the MGF stochastic network calculus performance bounds 
on backlog and delay were generally derived in (7) from 
Chernoff's theorem for arrivals and service random process 
that posses the MGFs M A (0, t) and M s (9, t) respectively. One 
of the strengths of stochastic NetCal lies in the characterization 
of end-to-end service by means of composition of the service 
provided at each hop, i.e. finding a single dynamic server char- 
acterization equivalent to multiple dynamic servers in series. 
For j] servers in tandem providing stochastically independent 
service, the MGF of the end-to-end system is bounded by [Th. 
1 from Q] for t > 



„(e,t)<(M Sl (0)*---*M Sri (0)) (t) 



(2) 



where * is the convolution operator known from system theory. 

A stochastic bound on the delay W such that P [W > d] < 
e assuming FIFO scheduling is derived in [7], with 
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B. MIMO with Spatial Multiplexing 

The M x 1 received signal vector y of a MIMO system 
with N transmit and M receive antennas is given as 




where s is a N x 1 transmitted signal vector, while n is a M x 1 
zero mean, complex Gaussian noise with independent and 
identically distributed real and imaginary parts and p is the av- 
erage SNR at the receiver. We further assume E[nn^] = A Im 
i.e. the noise corrupting the different receivers are independent 
where denotes the matrix conjugate transpose, N a is the 
noise power density and Im is the identity matrix of dimension 
(M). H is M x N MIMO channel matrix given as 
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where h mn represents the fading path gain from the nth 
transmit antenna to the mth receive antenna. We use the finite 
scatterers MIMO channel model (23), which is widely used for 



describing multipath environments, to describe the individual 
channel gains. Under this model, the channel path gain h m<n , 
assuming no line of sight (NLOS) component, is given as p4| 



a s A s exp 



8=1 



2tt 



(3) 



where iV a is the number of scatterers, a s is a uniformly 
distributed complex scattering coefficient for a scatterer s. 
D n s is the distance between the transmit antenna n and the 
scatterer s while D s , m is the distance between the scatterer s 
and the receive antenna m and A s is the attenuation constant 
of the communication path. 

The squared Frobenius norm of H, ||H||^ represents the 
total power gain of the channel and is given as 



M N 

EE 

m=l n— 1 



The meaning of the SNR p depends upon the normalization of 
the channel matrix [25 J. We apply the most commonly used 
normalization (2) 



H 



= N ■ M 



In this case, p is the SNR per receiver from all transmitters 
while p/N is the SNR per receiver per transmitter. 

The MIMO capacity under a given realization of channel 
matrix for uncorrelated equal power sources is given as (|2j 

P 



C = log 2 



det I 
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-HH 



(4) 



MIMO offers spatial multiplexing gain in which indepen- 
dent data flows are transmitted over different antennas to 
maximize the data rate. Spatial multiplexing takes advantage 
of the extra Degrees Of Freedom (DOF) made available due to 
the multiple antennas. A MIMO channel has up to min(N, M) 
DOF JT], Q. In this work, we quantify the spatial multiplexing 
gain by modeling the available DOF in a MIMO system using 
a discrete time Markov Chain. The available DOF, among 
others, depends upon the physical environment in which the 
MIMO system operates, which is represented by the channel 
matrix elements. We next formulate the problem for QoS 
analysis of MIMO system for the general case in which entries 
of H are formulated as in ||3}. 

III. The Analysis 

The discrete time block fading model of a MIMO system 
with N and M transmit and receive antennas respectively can 
be represented by a discrete time Markov chain. An immediate 
idea is to use 2 NxM states to describe the system. However, 
it introduces exponential explosion of the state space and 
consequently the computation complexity. 

In our analysis, we define states for the Markov chain based 
on the concept of degree of freedom (DOF). Particularly, for 
the N x M MIMO, K = min(N, M) + 1 states are adopted, 
where the ith state denotes that there are (i — 1) DOF available 
in the channel and has a stationary state distribution vector 7r. 




Fig. 1. MIMO with Spatial Multiplexing Mode 



The Markov chain evolves according to the states of the un- 
derlying spatial paths, which dictate the available DOF. These 
spatial paths form the sub-states of the Markov Chain. A par- 
ticular path h m n can be either {good} or {bad} (abbreviated g 
resp. b) under the GE model, depending upon its power which 
is in turn dictated by the number of scatterers N s in ([3]). The 
2x2 scenario including four spatial paths {tin, hi2, ^21, ^22} 
is depicted in Fig. [T] A combination {gggg} denotes that all 
the paths are in the good state. 

We assume that the path stays in the good state with 
probability l—p g b and in the bad state with probability l—pb g . 
The probability to go from the good state to the bad state is 
therefore p g b and from the bad state to the good state is pi, g . 
For a certain path the probability of being in the {bad} state, 
i.e. the block error probability for that path, is given as 



Pgb 



Pgb + Pb g 



(5) 



It is worth highlighting that for spatial multiplexing MIMO, 
whenever two paths targeting the same receiver antenna are 
{good}, such as {ggbb}, the receive antenna can only decode 
one of them. In addition, if more than one path with the 
same transmitter are {good}, such as {gbgb}, the transmission 
can be counted only once. Note that we are not analyzing 
automatic repeat request (ARQ) behavior. We assume that 
whenever all paths from a transmit antenna are {bad} the data 
blocks intended to that antenna are redirected to the other 
antennas. 

For the Markov chain whose state is based on DOF, the 
transmission rate in state i is given as 



r % = min{C Sti } 



(6) 



where C s .i is the average capacity of sub-state s of state i 
obtained by taking the average of C in Q. Here, can be 
considered as the minimum capacity achievable in the state i, 
which holds in the context of delay limited capacity | |26) . 

The rates r; are ordered in the diagonal matrix R(9) 
as diag(e 



9ri 



,e 



8r K 



). The first order capacity is given as 
X)i=i r i 7T i- Th e transition probability matrix is represented as 
Q, where its element pij denotes the transition probability of 
state i to state j. The MGF of the random process S described 
by such Markov chain is given for 9 > and t > by [8] 



M s (0,t) = 7r(R(-6»)Q)*~ 1 R(-6»)l 



(7) 




Fig. 2. States of the Markov Chain for N = 2 

where 1 is column vector of ones. 

Next we consider the simplest case of 2x2 MIMO operating 
in spatial multiplexing mode. A 2x2 MIMO has 2 degrees of 
freedom and is hence described by a three state Markov chain: 

. State 1: DOF= 0. 

. State 2: DOF= 1. 

. State 3: DOF= 2. 

Each state is an aggregation of the sub-states where the 
aggregation is performed on the basis of DOF. Fig. [2] depicts 
this aggregation. The rate supported by the state i with DOF 
i — 1 is the minimum rate achievable in this state. 

. State 1: bbbb. 

• State 2: ggbb, bbgg, gbgb, bgbg, gbbb, bgbb, bbgb, bbbg. 

. State 3: gggg, gggb, ggbg, gbgg, bggg, gbbg, bggb. 
The transition matrix Q of this three state process can be 
calculated using conditional probabilities and the law of total 
probability. The steady state vector follows accordingly using 
7r = 7r ■ Q, which is given here as 

7r lx3 = [w 4 4w 2 (l-w) 2 +4w 3 (l-w) 
(1 - uj) 4 + 4w(l - w) 3 + 2w 2 (l - oj) 2 } 

The corresponding rate matrix R is given as R(0) = 
diag(l, e 8r2 , e 3 ), where we assume that there is no workload 
processed in state 1. The method is easily extended to NxN 
MIMO systems. 

With |7]l decided, delay bounds for the MIMO system 
are readily obtained using stochastic network calculus results 
discussed in Sec III-AI 

IV. Numerical Results 

In this section we present numerical results for the NxN 
MIMO channels as described in the previous section for differ- 
ent A. The delay bounds presented are calculated numerically 

We use a parameterization 



based on the formulae in Sec II-A 



for the channel and for the traffic source to compute the 
numerical results for a NLOS scenario. 

We parameterize the MIMO channel according to 802.1 In 



the discrete time model of 31/is which is approximately the 
time needed to transmit a data block of 2312 Bytes over a 
600 Mbps channel, the maximum rate given in (27). As we 
consider scenarios with AxA antennas we obtain the capacity 
of each of the DOF states through an extensive simulation 



as described in section II-B We omit showing confidence 



intervals as these turn out to be very small. We set A s = 500 
to describe a good channel which employs that the NLOS 
component power is enough to support error free transmission 
while a bad channel is described by zero NLOS component 
power. We set A s — 1 Vs while the free space attenuation 
factor l/(Z?„ s _D s m ) is taken to be a constant such as in |24|. 
For the GE model, we set the transition probability from {bad} 
to {good} pi,g =0.1 and the probability from {good} to {bad} 
Pgb = 0.01 similar to fl6|. 

For the discrete time model we use a periodic source 
that generates arrival traffic. Such traffic source with pe- 
riod r produces a units of workload (data blocks) at times 
{Ut + nr, n = 0, 1, . . .} where U is the starting time which 
is uniformly distributed in the interval [0, 1]. The MGF of such 
a source is known (see e.g. (22)) as 



M 



,t) = e MAl (i 



using 40 MHz channels [27]. We choose a base time unit for 



for t > and 9 > 0. The framework described before 
enables the derivation of probabilistic performance bounds, 
such as for the delay, for a various number of traffic sources 
whenever the MGF exists. A collection of such MGFs for 
sources of different statistical properties can be found in p2~[ 
and references therein. In this work we focus on the periodic 
source model and exploit the enhancement in QoS, in terms 
of delay bound, delivered by the use of multiple antennas 
at the transmitter and receiver. The periodic source can be 
seen as a streaming source, generating a data blocks every 
time period r, which enables setting the arrival rate v. The 
source arrival rates are tuned to meet modulation and coding 
schemes (MCS) from the IEEE 802.1 In standard (27), i.e. rate 
240 Mbps corresponds to MCS- 13 for N — 2 antennas. We 
parameterize the source such that its period r is 10 time units 
and set the number of generated data blocks a to achieve a 
given arrival rate v. We fix e — 10~ 6 and v = 240 Mbps 
unless stated otherwise. 

Fig. [3] depicts numerical results on the stochastic delay 
bound as a function of the arrival traffic rate v for different 
violation probabilities e. Observe the degradation in terms 
of delay bound as the violation probability gets tighter. As 
discussed above an increase in the number of antennas leads 
in the best case to a linear increase in channel capacity. This 
increase capacity manifests itself in the decrease of the delay 
bound for fixed scenario conditions. This can be seen in Fig. |4] 
Further, observe the exponential tail decay of the delay bound. 
This is expected for the delay bound as it is derived using the 
MGFs and Chernoff's theorem. The figure shows additionally 
the change in the exponential decay factor, seen in the slope 
of the lines, which is introduced through the increase in the 
number of antennas N, respectively, capacity. 
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Fig. 3. Stochastic delay bounds for N = 2 
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Fig. 4. Stochastic delay bounds for different number of antennas N 

We next study the effect of fading speed on the delay bound. 
We start from the configuration of transition probabilities pb g 
and p g b given before and fix the corresponding block error 
probability ui for each spatial path as given by (j5J. We control 
the fading speed, under fixed w, by the mean time for the 
Markov chain of the path to change states twice. To this end 
we variate pb g such that a small value corresponds to slow 
fading and vice versa. 

Tab. [I] verifies the known notion that higher antenna con- 
figuration provide a higher capacity. Note that the ratio of the 
first order capacities of N = 3 to N = 2 stay unchanged over 
different fading speeds. 

The impact of the fading speed on the delay bound for 
different number of antennas (N = 2,3) is depicted in Fig. [5] 
It shows that the delay bound deteriorates with slow fading 

TABLE I 

First order capacity for different # of antennas 
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Fig. 5. Impact of fading speed on the delay bound 
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Fig. 6. Impact of SNR on the delay bound for N = 2 



significantly. This is due to the fact that a slow fading channel, 
which corresponds to small values of pb g , produces large burst 
errors due to its longer memory. The notable thing is that for 
slow fading the difference in delay bounds between N = 2 and 
N = 3 is large and for higher fading speeds, which corresponds 
to a channel nearly without memory, the profit from extra 
spatial paths diminishes. As the first order capacity given in 
Tab. [I] stays constant over the range of fading speeds, which is 
implied by the fixed to, this difference in delay bounds might 
be due to the higher order statistics of capacity. Furthermore, 
the impact of statistical multiplexing between the spatial paths, 
which is higher for N — 3, is apparent in the slow fading 
regime where the Markov chain changes states less frequently. 

Fig. [6] shows results on the relation of the delay bound and 
SNR. One can observe a rapid growth of the delay bound 
occurring for small SNR values. This can be explained by the 
fact that the SNR governs the capacity expression (|4]) such that 
for a constant arrival rate the decrease in SNR corresponds to 
an increase in utilization, which, in turn for high values, causes 
large delays to occur more often and stochastic delay bounds 
to grow quickly. 
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Fig. 7. End-to-end delay bounds for different number of antennas on each 
hop. 

For Fig. [TJ we consider a simple multihop scenario, where 
we assume statistical independence of the service provided by 
the individual links. This scenario describes an ad-hoc wireless 
network, where we use homogeneous links for simplicity, i.e. 
all links are operated by NxN MIMO in spatial multiplexing 
mode, here N € {2,3,4}. Note that the framework described 
in Sect. [II] enables the derivation of delay bounds for a 
heterogeneous setting, by using the corresponding MGFs. 
Fig. [JJ depicts the growth of the probabilistic delay bound in 
the number of wireless hops. For the MGF of the end-to-end 
service we use |2]). The first observation in Fig. [JJ is that the 
delay bounds scale at most linearly in the number of hops -q. 
This observation is in accordance with the scaling of O (77) 
derived for MGF NetCal first in |7j. The second observation 
is that the number of antennas TV changes the slope of the 
delay bound growth in the number of hops rj. The additional 
capacity, which is provided by the allocation of additional 
antennas in the MIMO system, causes the change in the slope 
of the bound scaling. 

V. Conclusions and Future Work 

In this paper we conduct a flow-level performance analysis 
of a MIMO system with spatial multiplexing. The analyti- 
cal framework we use is stochastic network calculus with 
moment generating functions. We model the block fading 
MIMO wireless channel as a Markov Chain and parameterize 
it to the application of IEEE 802.1 In. The states of the 
Markov chain are decided on the degree of freedom, which 
results in a significant state space reduction from a direct 
way of defining the states. We show the improvement related 
to spatial multiplexing on the flow-level probabilistic delay 
bound. We also present numerical results of the delay bound 
and quantify the impact of increasing the number of antennas. 
Furthermore, we show results for end-to-end delay bounds in 
the independent multi-hop wireless scenario and the impact of 
increasing the number of antennas characterizing the MIMO 
links. While in this work periodic source has been used in 
illustrating the results, the same methodology can be applied to 



any traffic source with known MGF. For multi-user scenarios, 
it is left as a future work. 
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